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Abstract 

We consider a mechanism of dark matter production in the course of first order 
phase transition. We assume that there is an asymmetry between X- and X-particles of 
dark sector. In particular, it may be related to the baryon asymmetry. We also assume 
that the phase transition is so strongly first order, that X-particles do not permeate 
into the new phase. In this the bubbles of old phase collapse, X-particles are 

packed into Q-balls with huge mass defect. These Q-balls compose the present dark 
matter. We find that the required present dark matter density is obtained for the 
energy scale of the theory in the ballpark of 1-10 TeV. As an example we consider a 
theory with effective potential of one-loop motivated form. 

1 Introduction. 

The idea that the baryon asymmetry and dark matter in the Universe may have common 
origin, is of considerable interest [1-20]. In particular, dark matter particles (X, X) may have 
their own conserved charge and be created, together with baryons, in asymmetric decays of 
heavier paricles [3-8, 15, 17,20]. In this scenario, the initial asymmetry in the dark matter 
particles is roughly of the order of the baryon asymmetry, 

nx-nx^ riB- (1) 

Assuming that X-X annihilation cross section is not particularly small, one obtains the 
estimate for the present dark matter mass density 



Px ~ mxTiB- 



(2) 



Thus, the correct value of px appears to require the mass of X-particles in a few GeV range. 

New physics, whose manifestation would be the existence of X-particles, may well be 
characterized by much higher energy scale, and the X-particle mass may grossly exceed a few 
GeV. One may wonder whether the co-production of the baryon asymmetry and dark matter 
can still work in this case. The estimate (2) shows that heavy X-particles are overproduced, 
so one needs a mechanism that makes the actual mass density of dark matter much lower 
than that given by (2). 

In this paper we consider a possible scenario of this sort. The idea is that X-particles 
may be packed into Q-balls. A Q-ball made of X-particles of total number Q typically has 
a mass ttiq, which is much smaller than Q ■ mx [21-25]. Hence, the mass density of the 
dark matter Q-balls is naturally well below the estimate (2). A mechanism that packs free 
particles into Q-balls applies to the Friedberg-Lee-Sirlin Q-balls [21-24] (as opposed to the 
Coleman Q-balls [25] explored in supersymmetric theories [26,27]) and is as follows [28]. 

Let us assume that X-particles obtain their mass due to the interaction with an addi- 
tional scalar field 0, so that mx = h(j), where h is the coupling constant. Let us also assume 
that there is the first order cosmological phase transition at some temperature Tc from the 
phase = to the phase = 0c 7^ 0, and, furthemore, that the X-particle mass is large in 
the new phase, /i0c ^ T^. Then X-particles get trapped in the remnants of the old phase, 
which eventually shrink to very small size and become Q-balls (see Fig. 1). 




Figure 1: Q-ball formation during the first-order phase transition from the phase = 
(white) to the phase 7^ (blue). 



We find that this mechanism indeed works in a certain range of couplings characterizing 
the model. Interestingly, the correct value of the present Q-ball mass density is obtained for 
X-particle mass in the ballpark of 1-10 TeV. This makes the scenario potentially testable in 
collider experiments. The Q-ball mass and charge are in the range mg ~ 10~^ — 10~^ g and 
Q ~ 10^^ — 10^^, respectively. In this respect our Q-ball dark matter is not very different 
from that discussed in the context of supersymmetric theories [26,27]. Phenomenology of 
our dark matter Q-balls is also similar to that of supersymmetric Q-balls [26,27], except that 
the latter may be stabilized by baryon number and hence eat up baryons. 



2 



The paper is organized as follows. We begin with a brief description of Q-balls (Section 
2). In Section 3 we discuss the creation of Q-balls in the course of the phase transition and 
relate the dark matter Q-ball parameters and their present mass density to the properties of 
the phase transition. We also consider the conditions of validity of our scenario. In Section 
4 we give a concrete example based on one-loop motivated form of the finite temperature 
effective potential and present the ranges of parameters for which our mechanism is viable. 
We conclude in Section 5. 



2 Q-balls. 

2.1 Q-ball configuration 

Q-balls are compact objects that exist in some models possessing a global symmetry and as- 
sociated conserved charge. One of the simplest models admitting Q-balls is of the Priedberg- 
Lee-Sirlin type [21-24]. Its Lagrangian is^ 

L = l{d,(l>r - U{<t^) + {d.xnd^x) - h'cPYx , (3a) 
C/(0) = A(0^ - , (3b) 

where is a real scalar field and x is a complex scalar field. The field % is meant to describe 
X-particles whose mass in vacuo equals 

mx — hv . 

These particles carry global charge, associated with the ?7(l)-symmetry x ~^ 6*'*X- The 
lowest energy state of large enough charge is a spherical Q-ball with (p — inside and (f) — v 
outside. At large Q, its size it! and energy E are determined by the balance of the energy of 
Q massless x-quanta confined in the potential well of radius R and the potential energy of 
the field in the interior, i.e., R and E are found by minimizing 

m) = ^ + Y^'uo , (4) 

where Uq — U{0) — U{v) — Xv^. Hence, the Q-ball parameters are 



"'^As it stands, this model has the discrete symmetry cf) — > —(f) and hence is not viable because of the 
domain wall problem. To get around this problem one can either assume that the discrete symmetry is 
explicitly broken or consider the field (p carrying some gauge or global charge. This qualification is irrelevant 
for what follows. 
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Note that the surface energy is proportional to Rq oc Q^^^ and therefore neghgible at large 
Q. The Q-ball is stable provided its energy is smaller than the rest energy of Q massive 
X-quanta in the vacuum (f) = v, 

rriQ < mxQ ■ (6) 
A Q-ball is a classical object, since its raduis is much larger than its Compton wavelength. 



2.2 Radius of cosmological Q-balls 

Assuming that Q-balls made of X-particles compose dark matter and that the X-asymmetry 
is related to the baryon asymmetry via Eq. (2), we can obtain an estimate for the present 
radius of a typical Q-ball already at this point. Indeed, we find from Eq. (5) that 

3mQ 

Now, since X-particles are packed into Q-balls, we have 

{nx - rix) _ uqQ 
s s 

where nq is the number density of Q-balls and s is the entropy density. Therefore 

^ = ^^^ = 3A^, (7) 
s s 

where Ab = ub/s = 0.9 x 10^^°, and we assume for definiteness that the X-particle asym- 
metry is equal to the quark asymmetry. Making use of the relation pdm = "m-Qng, we 
obtain 



to Pdm 

where the subscript " 0" refers to the present epoch. Hence, we obtain from Eq. (5) that the 
present radius is 

Rq = 47r^ ~ 6 • 10-1^ cm. (9) 

Pdm 

Thus, even though the Q-ball mass and charge depend on the parameters of the model, its 
typical radius does not. 



3 Q-ball production at the cosmological phase transi- 
tion. 

In this section we give a general description of the Q-ball formation without specifying the 
form of the finite temperature effective potential. We assume only that it has a minimum 
= at high temperature, and that the phase transition from = to 7^ is of the first 
order. 
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fc=^^. (12) 



3.1 Bubble nucleation rate 

Let Tc be the critical temperature at which the effective potential has two degenerate minima 
at (f) — and (p — (pc- Below this temperature, the new phase (j) — (pc has lower free energy 
density, 

VefMc) - K//(0) = -p < 0. (10) 

Thermal fluctuations lead to the creation of the bubbles of the new phase. In the thin wall 
approximation, the free energy of a bubble of radius R is 

4 

F{R) = --ttR^p + 47ri?V , (11) 

where a is surface free energy density. The extremum of (11) gives the free energy of the 
critical bubble, 

IGtt 

It, in turn, determines the bubble nucleation rate, 

„ r^A -Mil 

where k is a factor roughly of order 1. Let us introduce 

Tc-T 

^-^Tf^-^ (13) 

c 

and assume that this parameter is small. To the leading order in r, we have 

p = Pit , (14) 
(7 = (7o, (15) 

where pi and ao are constants independent of r. Thus 

r = «7;V^^, (16) 

where 



^=^5- (17) 



167r(7g 

As the temperature decreases, the bubble nucleation rate rapidly grows 
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3.2 Temperature of the phase transition 

To estimate the number density and charge of Q-balls produced during the phase transition, 
we need an estimate of the transition temperature. Since the duration of the phase transition 
is often much shorter than the Hubble time (this corresponds to r ^ 1), we neglect the 
cosmological expansion and write for the fraction of volume occupied by the old phase at 
time t [29] 

x{t) = exp[-A(t)] , 

where ^ 

A{t)= I V{t,t')r{t')dt' , 

Jtc 

V{t, t') = ^ [u{t — t')]^ is the volume, at time t, of a bubble of the new phase born at time 
t', and u is the velocity of the bubble wall. We make use of the standard relation between 
the Hubble parameter and temperature, H — T'^/Mpi, where Mpi — ^ and g'* hereafter 
denotes the effective number of the degrees of freedom at the phase transition temperature. 
We obtain for r <^ 1 

This integral is saturated near the upper limit of integration, and we get 

4 ^12 



~ kmM ^ 1 ^e"^ . 

The phase transition occurs when A is roughly of order 1, which happens in a narrow interval 
of temperatures around 

T. = AV2^-1/2 ^ (19) 



M%,\ T 



A 



where 



In 



, M* ^ ^ 



(20) 



with logarithmic accuracy. Our estimate is valid provided that r* <^ 1, i.e. 

A < L . (21) 

In what follows we assume that this is indeed the case, see also Section 4.2. 

Note that at the time of the phase transition, the bubble nucleation rate (16) is still 
small, 

(22) 



This is of course a consequence of the slow cosmological expansion. 
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3.3 Q-balls in the end of the phase transition 

We are now ready to estimate the volume from which X-particles are collected into a single 
Q-ball. This volume will determine the number denisty of Q-balls immediately after the 
phase transition and the typical Q-ball charge. One way to obtain the estimate is to notice 
that within a factor of order 1, this volume is the same as the volume of a remnant of the old 
phase in the midst of the phase transition, see Fig. 1. We estimate the size of a remnant 
by requiring that it shrinks to small size before a bubble of the new phase is created inside 
it. The lifetime of a remnant of size R is R/u, so the latter requirement gives 

Rim^r.!. 

Making use of Eq. (22) we obtain 

- ^P-- (23) 

Another way to estimate the volume that will shrink to one Q-ball is to consider bubbles of 
the new phase instead, and estimate the typical size of a bubble in the midst of the transition. 
At time t, the average bubble volume is 



R\t) = N-\t) f V{t,t')V{t')x{t') dt' , 

Jtn 



where N{t) = T{t')x{t') dt' estimates the number density of the bubbles. These integrals 
are evaluated in the same way as Eq. (18), and we get 

-RM^2n^-^ 

This gives the same estimate as (23), which demonstrates the consistency of the approach. 

From now on we use the following expression for the volume from which X-particles are 
collected into a single Q-ball: 

3 * ^ V 

where we inserted r* given by Eq. (19), and ^ is a parameter of order 1 which parmetrizes 
the uncertainty of our estimate. 

The number density of Q-balls immediately after the phase transition equals nQ{Tc) — 
V~^, and its ratio to the entropy density is 

nQin) ^ ^5 1_ 
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Note that for given values of couplings, this ratio is suppressed by (Tc/Mpi)^. Since one Q-ball 
contains all excess of X-particles in volume K, its charge is (again assuming X-asymmetry 
equal to quark asymmetry, cf. Eq. (7)) 

Q = nxV, = 3ABs(Te)K ■ (25) 

This charge is large, since it is proportional to {Mpi/T^)^. The X-particles are packed into 
Q-balls rather efficiently. 



3.4 Q-balls at present. 

Once the phase transition completes and Q-balls get formed, the ratio of their number density 
to entropy density stays constant and given by Eq. (24). With the Q-ball charge (25), its 
mass is found from Eq. (5), 



Hence, Q-balls are dark matter candidates provided that 



7.3.5V4.a^3/4j.3/4m;V4!LJl_^_, (26) 



rp3/4rrl/4 
"^QT^Q _ ^A3/4 -tc Up _ PdM 

s -^^^ M*,3/^ " So ' ^'^^ 

where pdm — 1 • 10~^ GeV • cm~^ and sq — 3000 cm~^ are the present mass density of dark 
matter and entropy density, respectively, and 

t-l/4 7-9/8 

/ = 19.0^ — -J. , (28) 

^ /l3/8^y%3/4' ^ ' 

is a combination of dimensionless parameters. As we see, the dependence on parameter ^ 
is weak, so we set ^ = 1 from now on. At this point we can make rough estimate for the 
relevant energy scale. Assuming Uq ~ T^, A ~ 1, ti ~ 0.1, gi* ~ 100 and L ~ 100, we get 
from Eq. (27) that should be in a ballpark of 

Te ~ 1 ^ 10 TeV. (29) 

As we pointed out in Introduction, this relatively low energy scale makes the scenario inter- 
esting from the viewpoint of collider experiments. 
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3.5 Validity of calculation. 

A particle physics model in which our mechanism can work must satisfy several requirements. 
One is the condition (21) or 

A « 41n . (30) 

Another is that X-particles do not penetrate into the new phase in the course of the phase 
transition. This is the case if their mass in the new phase is sufficiently larger than tem- 
perature. Quantitatively, we require that the mass density of remaining free X-particles is 
negligible compared to the mass density of dark matter Q-balls, 

{nx + nx) mx|(T=o) ^ Pdm 

S So 

Let us see what this condition means in terms of parameters. 

The dynamics of penetration of X-particlcs into the new phase depends on the bubble 
wall velocity and the strength of X-particlc intcrations with cosmic plasma. There arc two 
extreme cases: very slow motion of the bubble wall and very fast motion. Let us consider 
them in turn. 

Slow wall. 

Let the wall velocity be so small that there is complete thermal equilibrium for X-particles 
across the wall. Then the chemical potentials in the old and new phase are equal, and since 
X-particles in the old phase are massless, the chemical potential is negligibly small: indeed, 
in the old phase /iT^ ~ {nx — n-x) ~ A^T^, hence fi/T ~ A^. The number density of 
X-particles in the new phase is given by the equilibrium formula for non-relativistic species, 
and we have 

^^ + ^^ ^ ,-1 . 2 f ^Y\-m./T ^ 1 f!!^y/\-m./T (33) 

s V 27r y \2t:T ) ' ^ ' 

where all quantities, including mx, are evaluated at Tc- 
Fcist wall. 

In the opposite case of fast wall, all X-particles that penetrate the new phase stay there. 
Flux of X-particles and X-particles, whose momentum normal to the wall exceeds mx is 



2 

jxx — 



J J mx 



This integral is saturated at pr near mx and pi <C pt- Hence, we write 

Pt+Pl^Pt+ 7^ — = Pt + 



2pT 2mx 
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and obtain 

Jxx = -^mxT^e--/^. 

As the wall moves, its radius increases by dR — udt, and the number of pernetrated particles 
is 

jxxSdt = JxxW'^ATrR'^dR. 
So, the number density in the new phase is in the end 

nx + nx = u'^ixx- 

This gives the estimate similar to that in the slow wall case, except that instead of particle 
velocity yjT jmx it involves wall velocity. 

If XX-annihilation is switched off in the new phase, then the ratio [nx + nx)/s stays 
constant after the phase transition. For any of the above cases, and hence for all intermideate 
ones, up to logarithmic corrections we get from Eq. (31) 



mx 



> In f (33) 



or 



This condition does not apply if XX-annihilation is efficient in the new phase. If Eq. (34) 
does not hold, the behavior of X-particles in the new phase is similar to that of WlMPs (note 
that the X-particle abundance in the new phase just after the phase transition either coincides 
with, or exceeds the equilibrium abundance). In that case the condition (31) implies that 
the XX-annihilation cross section exceeds the standard WIMP annihilation cross section. 
Since the energy scale inherent in the model is high (see Eq. (29)) the latter scenario is not 
particularly plausible. In what follows we assume that the inequality (34) holds. 

4 Example. One-loop motivated effective potential. 
4.1 Critical temperature and Q-ball parameters. 

As an example, let us consider effective potential of a particular, one-loop motivated form, 

[/((/), T) = a{T^ - T^,)<l>' - ^T4>^ + A0^ (35) 

where Tc2, a and 7 are parameters depending on particle physics at temperature T, and we 
neglect temperature corrections to the quartic self-coupling. We treat Tc2-i ot-i 7 and A as 
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parameters of the model, but keep in mind the relation q;T^2 — SAv^, which follows from Eq. 
(3b). 

We assume in what follows that 97^ < 32q;A. Then at high temperatures, the effective 
potential has only one local minimum at = 0. As the Universe cools down, the minimum 
at 7^ develops. It becomes deeper than the minimum at = at T < Tc, where 

In what follows we need the expression for the height of the scalar potential and mass of the 
0-particle, both at zero temperature, in terms of the parameters entering (35) and (36), 

^_ (4aA-7T 4 /4aA-fy/% 

^ 64A3 ' ^^=i^^j ^'^^ 

In this way, we trade the parameter v in the original Lagrangian (3a) for the parameters 
relevant for the phase transition. 

Let us now rewrite the results of Section 3 in terms of the parameters used in (35). At 
small T = one has 

J c 

^- ^^(4aA-7!)^,^^ (38) 



8A3 

„3 



Thus 



and 



a^l^m^W-^Tf. (39) 

81V2A3/2(4q;A-72)' ^ ^ 

* V T^L^'^ ) V9 ■ 2^/^ T2 A3/4(4aA - 7^) ) ' ^ ' 

where L is defined in Eq. (20). We find from Eq. (27) the present Q-ball dark matter mass 
density 

POM = XpAf .oM^r'/ V'/'^^-'/'t-^'/'A-^/^^ (4aA - ^^f" T'J'L'/', (42) 

where 

Kp = 215/1637/451/4^1/8 ^ 22.5. 

Equating it to the actual dark matter density we obtain the critical temperature in terms of 
other parameters, 

n = i^Tpi/L A//^.-'/'M;,^/^^y7^3/7^15/14^3/28 (4^^ _ ^2) -5/7 ^ ^_9/,4^ (43) 
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where 

ii'y = 2-15/283-15-1/7^-1/14 ^Q_^7^ 

In further computations we use Tc = 10 TeV in the argument of logarithm (20), see Eq. 
(29). Finally, Eqs. (25) and (26) give for the Q-ball parameters 

Q = KQp-T<'4"^'M*^,''/'g^'/' 7=^°/' u''/' c?'l'\'l' (4aA - 7^) ~' L-'^l\ (44) 
m<, = i^.,pB^i^A^/^4/^M; -/^,//^z.-/V°/^A-«/^ (4aA - ^^y'" L-^^l\ (45) 

where 

Kq = 2^9/7^26/73-45-4/7 ^ Q Qg3^ 

K^^ = 2i=^/^3-55-^/V6/^ ^ 0.42. 
The present number density of Q-balls is of course equal to Pdm/^q- 

4.2 Parameter space. 

In Section 3.5 we pointed out two conditions that the model should obey. In terms of the 
parameters of the effective potential the condition (30) takes the following form (see Eq.(40)) 

«41nf^Ul20, (46) 



We find 



81v^A3/2(4q;A-72)2 V Tc 

The concrete form of the condition (34) is obtained by noticing that 

<t>c= mT=n^^n. (47) 

> 25. (48) 

We also assume that the quartic self-coupling of the field is not particularly small, and 
impose a mild constraint motivated by naturalness argument. 

The conditions (46), (48) and (49) are actually quite restrictive. In particular, they require 
that the 0-x coupling is rather strong. For /i ~ 5, the available region in the parameter 
space is fairly large, as shown in Fig. 2. This region becomes considerably smaller already 
for /i ~ 3, see Fig. 3. 
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Figure 2: The region in the parameter space consistent with the constraints (46), (48), (49) 
for h = 5. Dots are the points for which the mass, critical temperature and Q-ball parameters 
are hsted in Tables 1 and 2. 



We scanned the available regions and found the following range of masses and critical 
temperatures, at which our mechanism does work 

m,^ G 3 ■ 10^ ^ 1 ■ 10^ GeV, 

Te G 1 ■ 10^ ^ 3 ■ 10^ GeV. 
For the Q-ball parameters we obtained the range 

mg G 3 ■ 10"^ ^ 3 • 10"^ gram, 

Q G 10^^ ^ 102^ 

nqel- 10"^^ ^ 3 ■ 10~^^ cm'l 

Three particular examples, corresponding to the points in Fig. 2, are listed in Tables 1 and 
2 for two values of the bubble wall velocity. 
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Figure 3: Same as in Fig. 2, but for h = 3. 



Table 1. Q-ball parameters for particular values of couplings (three left columns) and 
bubble wall velocity u = 0.03. 









m<^ ,GeV 


Te,GeV 


Q 


mQ,gram 


uq, cm"^ 


a = 5 


7 = 0.4 


A = 0.04 


1.9 X 10^ 


4.8 X 10^ 


8.7 X 10^9 


4.2 X 10-^ 


4.5 X 10-26 


a = 0.5 


7 = 0.02 


A = 0.002 


8.1 X 10^ 


6.1 X 10^ 


7.2 X 10^9 


2.3 X 10-^ 


8.1 X 10-26 


a = 0.5 


7 = 0.004 


A = 0.0004 


3.7 X 10^ 


2.6 X 10^ 


1.9 X 10^9 


5.6 X 10-^ 


3.4 X 10-25 



Table 2. Same but for u = 0.3. 









,GeV 


Te,GeV 


Q 


mQ,gram 


nQ,cm-3 


a = 5 


7 = 0.4 


A = 0.04 


4.9 X 10^ 


1.2 X 10^ 


3.9 X 1021 


1.9 X 10-3 


1.0 X 10-27 


a = 0.5 


7 = 0.02 


A = 0.002 


2.1 X 10^ 


1.6 X 10^ 


3.2 X 1021 


1.0 X 10-3 


1.9 X 10-27 


a = 0.5 


7 = 0.004 


A = 0.0004 


9.4 X 10^ 


6.7 X 10^ 


8.2 X 1020 


2.5 X 10"^ 


7.8 X 10-27 



5 Conclusion. 



We considered a mechanism producing Q-balls in the course of the first order phase transition 
and described it quantitatively. As we have seen, this mechanism efficiently packs massive 
stable particles, which otherwise would be overproduced, and thus drastically reduces the 
mass density. 
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Using a well-studied model of Friedberg-Lee-Sirlin Q-balls, we obtained formulae for the 
dark matter properties, such as charge, mass and concentration of dark matter Q-balls, as 
well as masses of scalar fields, depending on the parameters of the transition. As an example, 
we considered this mechanism in a theory with the effective potential of one-loop motivated 
form. 

We have seen that the main parameter is the temperature of the phase transition, whose 
adjustment yields the right value of the present dark matter density. A remarkable property 
of the mechanism is that for a wide range of parameters, independently of particle physics 
and effective potential model, the estimate for the energy scale is Tc ~ 1-10 TeV. 

The main requirement for effective packing of particles into Q-balls is that the phase 
transition is strongly first order. For one-loop motivated effective potential this implies 
strong (j)-x coupling and gives the main constraint on available parameter space. However, 
for other models this is not necessaraly the case, since there are other ways to make the first 
order phase transition strong enough (see, e.g., [34,35] and references therein). 

Our study was motivated by the possibility that X-particle asymmetry is related to the 
baryon asymmetry. However, this is optional. The mechanism described works in the same 
way, with 3Ab replaced by Ax, if rix is considered as yet another free parameter. 
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